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DECOMPOSITION  AND  FICTITIOUS  DOMAINS  METHODS 
FOR  ELLIPTIC  BOUNDARY  VALUE  PROBLEMS 

S.  V.  Nepomnyaschikh' 


Abstract. 

Boundary  vzdue  problems  for  elliptic  second  order  equations  in  three-dimensional 
domains  with  piecewise  smooth  boundaries  are  considered.  Discretization  of  the 
problem  is  performed  using  a  conventional  version  of  the  finite  element  method  with 
piecewise  linear  basis  functions.  The  main  purpose  of  the  paper  is  the  construction 
of  a  preconditioning  operator  for  the  resulting  system  of  grid  equations.  The  method 
is  based  on  two  approaches:  decomposition  of  the  domain  into  subdomains  and  using 
a  new  version  of  the  method  of  fictitious  domains.  The  rate  of  convergence  of  the 
corresponding  conjugate  gradient  method  is  independent  of  both  the  grid  size  and 
the  number  of  subdomains. 


^  This  work  was  supported  by  the  Sloan  Soviet  Visitors  Program  grant  number 
89-4-8. 


1.  Introduction.  In  this  paper,  we  consider  preconditioning  operators  for  the 
system  of  grid  equations  approximating  the  following  boundary  value  problem: 

-  E?,j=i  g~^'^^^^'9~  +  ao(x)u  =  /(x)  ,      x  G  Q  , 

u{x)  =  0  ,  X  e  To  ,  (1-1) 

du 

—  +  a{x)u  =  0  ,      X  G  Fi  . 

uu 
Here,   — —  is  the  derivative  in  the  conormaJ  direction.     We  assume  that  fi  is  a 

oN 
polyhedron.  Let  Q  be  a  union  of  n  nonoverlapping  subdomains  Q,, 

n 

n  =  u  n, ,   Q.  n  fij  =  0 ,  t  ^  i 
1=1 

where  Cl,  axe  polyhedrons  with  diameters  on  the  order  of  H.  Let  Q^ 

1=1 

be  a  regular  triangulation  of  Q.  which  is  characterized  by  a  parameter  h. 

We  use  H^{Q,,ro)  to  denote  a  subspace  of  the  Sobolev  spaces  H^{Q)  of  real 
functions 

H\Q,  To)  =  {ve  H\n)  I  v{x)  =  0  ,  X  G  To} 

and  we  introduce  the  bilinear  form 

a{u,v)  =   /  (  y^  a,,(x)— — — \- aQ{x)uv)  dx -\-  /    a{x)uv  dx 

Jn     ~^^  oxj  ox,  Jr^ 


and  the  linear  functioned 


^{v)  =   f  f{x)vdx  . 


We  assume  that  the  coefficients  of  the  problem  (1.1)  are  such  that  a{u,v)  is  a 
symmetric,  coercive  and  a  continuous  form  in  H^{Q,,Tq)  x  i?^(f),ro),  i.e. 

a(u,u)  =  a(t;,u)  ,  Vu,t»  G  H^(^,Tq)  , 

ajid  that  the  linear  functional  i{v)  is  continuous  in  H^{Q,Tq)^  i.e., 

\t{v)\<^i\\v\W(u)  . 

Here  /xq,  /ii,  //  are  positive  constants.  A  weak  formulation  of  (1.1)  is: 
Find  u  G  H^{Q,,To)  such  that 

a{u,v)  =  e{v),      yveH\^,ro)  ■  (i-2) 


Denote  by  W  the  space  of  real-valued  continuous  piecewise  linear  functions. 
Using  the  finite  element  method  we  can  pass  from  (1.2)  to  the  linear  algebraic 
system 

Au  =  f  ■  (1.3) 

Our  purpose  is  the  construction  of  a  preconditioner  B  for  the  problem  (1.3)  such 
that  the  following  inequalities  are  valid: 

ci(Bu,u)  <(.4u,u)  <C2(Bu,u)  ,  \/ueR^.  (1.4) 

Here  A''  is  the  dimension  of  W,  the  positive  constants  ci,C2  axe  independent  of  h 
and  H,  and  the  multiplication  of  B~^  by  the  vector  can  be  at  a  low  cost. 

As  a  rule,  the  most  efficient  preconditioners  for  solving  boundary  value  problems 
in  domains  with  complex  geometry  can  be  constructed  by  simplifying  geometry  of 
the  original  domain.  We  use  two  approaches:  decomposition  of  the  domain  into 
subdomains  [5-7]  and  a  version  of  the  method  of  fictitious  domains;  cf.  [8-10]. 

The  remainder  of  the  paper  is  organized  as  follows.  In  section  2,  we  describe  ab- 
stract results  which  axe  useful  in  constructing  the  preconditioning  operator.  Using 
Lemma  2.1  we  spht  the  original  space  of  grid  functions  into  a  sum  of  subspaces  and 
using  Lemma  2.2  and  Lemma  2.3,  on  the  fictitious  domains  method,  we  give  equiv- 
alent norms  in  subspaces.  In  section  3,  we  construct  a  decomposition  of  the  original 
mesh  space  into  a  sum  of  two  subspaces.  The  functions  of  the  first  subspace  axe 
equal  to  zero  at  the  boundaries  of  subdomains.  The  second  subspace  corresponds 
to  the  values  of  mesh  functions  at  the  boundaries  of  subdomains  extended  inside 
with  conservation  of  the  norm.  The  second  subspace  has  a  complex  structure,  and 
to  simplify  it,  we  decompose  it  into  a  sum  of  a  coarse  space  and  local  subspaces  [1- 
3,9,11].  In  section  4,  using  the  abstract  results,  we  design  preconditioning  operators 
for  these  subspaces.  This  design  uses  the  fictitious  domain  (space)  method  to  over- 
come two  obstacles:  complex  geometry  of  the  substructures  and  chaotic  distribution 
of  the  nodes  of  the  triangulation. 

2.  General  Framework.  In  this  section,  we  outline  an  abstract  method  for 
constructing  a  preconditioning  operator  by  spHtting  the  original  Hilbert  space  into 
a  sum  of  subspaces.  Here,  we  completely  follow  [5-10]. 

Lemma  2.1.  Lei  the  Eilbert  space  H  with  the  scalar  product  {u,v)  be  split  into 
a  sum  of  subspaces 

H  =  H,+H2  +  ---  +  Hm  , 
let  A:  H  —*  H  be  a  linear,  self -adjoint,   continuous  and  positive  definite  operator, 
and  P,  :  H  —*  H,,  i  =  1,2, .. .  ,m,  be  the  orthogonal  projections  with  respect  to  the 
scalar  product  (u,f)^  generated  by  the  operator  A 

(u,  v)a  =  {Au,v)  . 

Assume  that  positive  constants  a  and  0  exist  such  that  for  any  element  u  ^  H  there 
exists  u,  G  H,  such  that 

Ui   +  U'2  +  ■  ■  •  +  Um   =  U 


Oi{  ("1,  "l)>l  +  ("2,  U2)a  H f-  {Um,Um)A)   <  (",  ")/l 

and 

(  (Pi   +  P2  +  •  •  •  +  Pm)u,  u)a   <  0{U,  u)a 

for  any  u    G    H.      Assume  further  that   there   are   selfadjoint   operators  S,,    i    = 
1,2, ...  ,m,  such  that 

Bi  :  H  —^  Hi  , 

ci{BiU,u)  <  (Au,u)  <  C2{BiU,u)  ,     Vu  G  Hi  . 

Then, 

aci(A-^u,u)  <  (B-^u,u)  <  ^C2(>l"^u,u)  ,      WuG  H  , 

where  B~^  =  B^  +  B2  +  •  •  •  +  5|J^,  an^  Bf  is  a  pseudomverse  of  B,. 

Using  Lemma  2.1,  we  will  defiBe  the  preconditioner  in  the  following  way.  Split 
the  original  space  of  grid  functions  W  into  a  sum  of  subspaces,  each  having  a 
definite  structure  which  maJces  it  possible  to  determine  easily  invertible  equivalent 
norms  (i.e.,  operators  B,)  in  these  subspaces.  The  following  lemmas  will  be  used  in 
constructing  operators  J5,,  which  give  equivsdent  norms  in  subspaces,  and  for  which 
the  multiplication  of  B^  by  the  vector  can  be  performed  efficiently. 

Lemma  2.2.  LetR"^  and  R"^  be  real  Euclidean  spaces,  and  E  and  S  be  symmetric 
positive  definite  matrices  of  order  m  and  n,  respectively.  Denote  by 

(<^,V>)s    =(E^,r/.), 
{u,v)s    ={Su,v) 

the  scalar  products  generated  by  these  matrices  in  i?*"  and  i?",  and  let  T:  R^  — +  i2", 
be  a  linear  operator  such  that 

for  any  4>  G  R"^-  Here,  a  and  0  are  positive  constants.   Let 

C  =  TE-'T'  , 

where  T*  is  an  operator  adjoint  to  T  with  respect  to  the  Euclidean  scalar  products 
in  R""  and  iJ".   Then, 

a{C'^u,u)  <  (u,u)s  ^  0{C'^u,u) 

for  any  u  £  E  =   Im  T.  C"*"  is  a  pseudo-inverse  of  C. 

Lemma  2.3.  Let  Ho  and  H  be  Hilbert  spaces  with  the  scalar  products  (uo,t'o)//o 
and  (u,t')//,  respectively,  and  let  A  and  B  be  selfadjoint  positive  definite  and  con- 
tinuous operators  in  the  spaces  Hq  and  H: 

A:  Ho-^  Ho  ,  B:  H  -^  H  . 


Let  R  be  a  linear  operator  such  that 

R:  H  — *  Hq  , 

{ARv,Rv)ho  <cr{Bv,v)„ 
for  any  element  v  E  H .  Furthermore,  let  there  exist  an  operator  T  such  that 

T:  Ho  —^  H  ,  RTuq  =  uq  , 

ct{BTuq,Tuq)h  <  (Auo,Uo)//o 
for  any  element  Uq  G  Hq.  Here,  cr  and  ct  are  positive  constants.    Then 

ct{A~'^uo,uo)  <  {RB~^R'uo,uo)ho  <  cr(A"'uo,uo) 

for  any  element  uq  G  Hq.  Here,  R'  is  the  operator  adjoint  to  R  with  respect  to  the 
scalar  products  (uo,Uo)//o  '^'^'^  ("i^)//  such  that 

R':  Hq  — >  H  , 

{R'uq,v)h  =  iuo,Rov)Ho  ■ 

3.   Decomposition  into  Subspaces.  To  design  the  preconditioning  operator 
B,  we  will  completely  follow  [9]  and  decompose  the  space  W  into  a  sum  of  subspaces 

W  =  Wo  +  W^  . 

To  this  end,  divide  the  nodes  of  the  triangulation  Q.^  into  two  groups:  those  which  he 
inside  of  Q,''  and  those  which  lie  on  boundaries  of  fif .  The  subspace  Wq  corresponds 
to  the  first  set.  Let 

1=1 

Wo    =  {u^  eW\u^{x)  =  Q  ,xeS^}  , 
Wo..     =  {tz'' e  Wo|u''(x)  =  0  ,  xefi,'*}  ,      z  =  l,2,...,n  . 

It  is  clear  that  Wq  is  the  direct  sum  of  the  orthogonal  subspaces  Wq,,  with  respect 
to  the  scalar  product  in  H^{Q.): 

Wo  =  Wo,x®---®  Wo,„  . 

The  subspace  W\  corresponds  to  the  second  group  of  nodes  Q''  and  can  be  defined 
in  the  following  way.  First,  define  W1/2  which  is  the  space  of  traces  of  functions 
from  W  on  S^: 

To  define  the  subspace  W\ ,  we  need  a  norm-preserving  extension  operator  of  func- 
tions given  at  S^  into  Q.^ .  Let  us  consider  the  subdomain  Qj*.  Let  {sx,S2,n)  be  a 
near-boundary  coordinate  system  [7,8]  which  is  defined  in  a  ^-neighborhood  of  d^^. 


Here  5i,  ^2  define  a  point  P  at  SQf  and  n  is  the  distance  between  the  given  point 
and  5Qf .  Set 

u(^i,^2,n)  =  (l- j) hrB—,[  4>{x)dx,  (3.1) 

0    Tcvezs  K{F,n)  JK(P,n) 

P  e  aJ2f  ,  0  <  n  <  6  , 

where 

A'(P,n)  =  {(5;,S2,0)e5Qf  ||5;-5.|  <n  ,  z  =  1,2}  . 

The  function  u  is  extended  by  zero  in  the  rest  of  Q,.  Using  Hardy's  inequality  [4], 
it  is  not  difficult  to  see  that  Ms  a  norm-preserving  operator,  i.e. 

Il"ll//'(fi.)  =  IK<?!'ll//'(n.)  -  (^A\<t>\\HU^^u,)  ,  V^  €    , 

where  Ci  is  a  positive  constcint.  Using  the  auxihaxy  mesh,  which  is  topologically 
equivalent  to  a  uniform  mesh,  we  can  define  the  finite-element  analogy  th  of  the 
operator  t  from  (3.1),  such  that 

\\tK<t>*'\\HHU,)  <  C2U''\\HU^9n,)  ,       ^4>'  (3-2) 

The  cost  of  the  multiplication  of  t  and  t'  by  a  vector  is  proportional  to  the  number 
of  mesh  domain  nodes.  Here  C2  independent  of  h,  H  and  t'  is  an  adjoint  operator 
to  t  with  respect  to  the  Euclidean  sccJar  product  (see  [8,9]  for  more  details).  Set 

It  is  obvious  that  W  =  Wq  +  Wi  and  according  to  (3.2)  this  decomposition  of 
W  satisfies  Lemma  2.1  with  a  independent  of  h  and  H . 

According  to  Lemma  2.2,  to  define  the  preconditioning  operator  Bi  for  the 
subspace  Wi  we  need  to  define  the  preconditioning  operator  P1/2  which  generates 
the  norm  in  the  space  Wi/2-  The  set  5''  has  a  complex  structure.  Therefore  we  will 
split  W-yi2  into  subspaces  associated  with  the  following  substructures.  Let  S^  be  a 
union  of  k  nonoverlapping  triangles  T,: 

h  _ 
5''  =  U  T,  .      T,  n  T,  =  0  ,  ?■  /  i  .  (3.3) 

«=i 

We  assume  that  k  and  n  are  of  the  same  order  of  magnitude.  Split  the  space  W112 
into  a  sum  of  subspaces 

W„2  =  Vo  +  V'l  +  •  •  •  +  14  ,  (3.4) 


where  Vq  is  the  coarse  space  which  consists  of  continuous  functions  hnear  on  the 
triangles  T,,  i  =  1,2,...,/:,  and  V,,  i  =  1, . . . ,  A;  correspond  to  T,  and  are  defined 
below.  Denote  by 

^1/2..  =  {<^?  €  w^i/2, 1  <^K^)  =  0 ,  X  G  5r.} . 

For  any  triangle  Ti,  we  define  the  explicit  extension  operator  r^: 

r,:  1^1/2,  -.  ^1/2  • 
Denote  by  T/  the  union  of  adjacent  triangles: 

i;'  =  _y    T,. 

Let  (f)^  be  a  function  given  on  T,.  We  will  define  the  extension  of  the  function  <f)^ 
using  two  tricks:  reflection  and  semi-reflection.  First  let  us  define  reflection  <p  of 
the  function  (p^.  Let  T,  be  a  triangle  from  7^',  which  has  a  common  side  with  T, 
(see  Fig.  3.1a). 


Figure  3.1 

Let  u,  V,  w,  u  +  V  +  w  =  1  he  the  barycentric  coordinates  on  the  triangle  T,. 
Then  any  point  x  from  T,  can  be  represented  in  the  following  form 

X  =  u-  Pi  +v  ■  P2  +  W  ■  P3 

where  Pi,P2,P3  are  the  vertices  of  the  triangle  T,.  Set 

4>(y)  =  4>{u-Pi  +  vP2  +  w-P^)  =  <t>{c)  =  <f>{u-Pi+vP2  +  w  P3) 

Now  let  us  define  semi-reflection  (f)  of  the  function  <f>  given  on  T,  for  the  case  when 
a  tricingle  Tj  has  a  vertex  common  with  the  triangle  T,  (see  Fig.  3.2). 


Fig.  3.2 


Set 

4>{z)     =4>iu-Pi  +  v-P^  +  w0.5-  (Pj'  +  ^3) ) 

=  4>{y)  =  4>{u-Pi+v-P^  +  w  0.5  •  (P2'  +  P^) ) 

=  4>{x)  =  <f>{uPi+vP2  +  w-0.d-  (P2  +  P3) )  . 

It  is  clear  that  we  can  define  the  "semi-reflection"  for  the  case  when  a  triangle  T, 
has  a  side  common  with  the  triangle  Ti  too.  Combining  both  approaches,  we  can 
define  an  extension  ^^  of  the  fimction  <^''  onto  T,.  Using  the  piecewise  linear  cutoff 
function  ^,  defined  by 

'  1  ,        xeT, 

0 ,     x^t:  , 

we  can  define  the  extension  operator  r,  by 


Ux)  = 


r,r  =  i,^ 


Set 


V,  =  T,Wi/2,.  +  W„2,i  ,      1  =  1,2,...,/:.  (3.5) 

Using  Lemma  4.1  from  [9]  it  is  easy  to  see  that  the  decomposition  (3.4),   (3.5) 
satisfies  Lemma  2.1  with  a  independent  of  h  and  H. 

4.  The  Fictitious  Space  Method.  In  this  section  we  will  design  a  precon- 
ditioning operators  for  the  subspaces  Wq,,,  i  =  1,2,  ...,n,  and  K,  i  =  0,1,...,  A; 
constructed  in  Section  3  using  Lemmas  2.2  and  2.3. 

There  axe  two  obstacles  to  using  the  finite-element  approximation  of  the  Laplace 
operator  as  a  preconditioning  operator  for  Wo,t-  The  first  is  related  to  the  geometry 
of  fi,  while  the  second  relates  to  the  chaotic  distribution  of  nodes  of  the  triangulation 

Q^ 

Let  the  domain  Q,  be  embedded  in  the  cube  E,.  We  extend  the  triangulation 
Q,'^  to  the  entire  cube  U'-  in  such  a  way  that  H-*  is  a  quasi-uniform.  We  will  denote 
by  I'r(nf )  the  space  of  continuous  functions  which  are  linear  on  each  simplex  of  the 


triangulation  n|',  and  by  An*-  the  approximation  of  a  symmetric  elliptic  operator 
in  VK(n^).  According  to  Lemma  2.3,  we  need  a  continuous  operator 

Ru,:W{U':)-^Wo,  . 

Let 

Rn,  =  In,  —  ih  ■  hn,  , 

where  In,,  /an.  are  the  trace  operators  from  11,''  onto  Q,'',  dO.'^  ,  respectively.  The 
operator  th  is  the  extension  operator  which  was  described  in  section  3.  From  the 
continuity  of  /n, ,  /an,  and  </,,  we  obtain  the  continuity  of  R\\,-  We  can  define  the 
operator 

as  an  extension  by  zero.  It  is  easy  to  see  that  the  hypotheses  of  Lemma  2.3  are 
valid.  Hence  there  exist  positive  constants  c\  and  C2,  independent  of  h  and  H  ^  such 
that 

for  any  u''  G  Wo,,. 

Further,  along  with  the  mesh  H^  in  the  cube  H,  we  will  also  consider  an  auxiliary 
uniform  mesh  Q'^  with  a  mesh  step  ho',  czh  <  ho  <  /imin/v3.  Here  hmi„  is  the 
length  of  the  minimal  edge  of  simplices  of  the  triangulation  Hf,  and  constant  c^  is 
independent  of  h  and  H.  Denote  the  nodes  of  the  mesh  Qf  by  Zj^k.i  =  (-Yj,  Yk,  Z(), 
j,  /:,  £  =  0, 1, . . . ,  M,  and  denote  the  cells  of  the  mesh  Q'^  by 

Qe.j,k  =  {(a:,  y,  z)\X,<x<  X,+i  ,Yk<y<Yk+i,Ze<z<  Z,+i}  , 
£,j,it  =  0,l,...,M. -1  . 

On  the  mesh  Qj',  we  will  consider  the  space  W{Q'l)  of  mesh  functions  U{Zj_k,()-  Let 
us  define  the  operator  Rq,: 

RQ,:W{Q^)-.WiIl^) 

which  introduce  a  correspondence  of  each  function  U{Z-,^k,e)  <  ^{Q'^)  to  a  function 
u^  6  W(nf )  in  the  following  way.  Let  Zm  be  a  node  of  the  triangulation  Hf  and  let 
Zm  e  Qj,k.(-  Let 

u'izrn)  =  U{Z,,k,()  . 

Note  that  by  the  assumptions  on  ho  at  most  one  node  z^  of  the  triangulation  H, 
can  belong  to  the  cell  Qj,k,t-  The  operator  Tq^, 

Tq,:W{^'1)-^W{Q'^) 

is  now  defined  as  follows.  If  the  cell  Qj,k,t  contains  a  node  2„,  of  the  triangulation 
nf,  we  set 


10 


At  the  other  nodes  of  the  mesh  Qf  the  function  U{zj^k,()  can  be  defined  in  a  relatively 
arbitrary  way,  for  instance,  as  follows.  Let  the  node  Zj^k,e  belong  to  the  simplex  r^ 
of  the  triangtilation  11^  with  the  vertices  z^^ ,  •  •  • ,  ^mi  •  Set 

U{z,,,,e)  =  \in'{zm,)  +  ---  +  u\zr,,),). 

It  can  be  easily  shown  that  the  operators  Rq,,  Tq,  defined  in  this  way  satisfy  the 
hypotheses  of  Lemma  2.3  and  that  the  corresponding  constants  do  not  depend  on 
h  and  H.  Thus,  there  exist  positive  constants  C4,  C5,  independent  of  h  and  H,  such 
that 

C4||u''||^Mn.)<((^n.Hg.(-Ag.)-'i?Q.i2nJ-'u,u)<C5||u''|l^:(„.) 

for  any  u^  G  Wq,,.  Here  Aq,  is  the  mesh  approximation  of  a  symmetric  elliptic 
operator  in  W{Q'^).  The  preconditioning  operator  Bq,  for  Wq,,  can  now  be  defined 
as  an  extension  by  zero  of  the  operator  Ru,Rq,{—^q,)~^ Rq^R^^-  Set 

We  now  define  the  preconditioning  operator  B^  for  the  subspace  Wi .  According  to 
Lemma  2.2,  we  can  define 

Bt  =  tH^-'t', 

where  the  operator  E  generates  an  equivalent  norm  in  the  trace  space  W1/2  and  the 
extension  operator  t^  was  described  in  section  3.  To  define  the  operator  S  we  will 
use  the  decomposition  (3.3).  Let  us  first  consider  the  coarse  subspace  Vq.  Let  Wq  be 
a  mesh  space  which  corresponds  to  the  values  of  the  functions  <f>  G  Vq  at  the  nodes 
of  the  triangles  T,  from  (3.3).  Define,  in  the  space  Wq,  the  symmetric  operator  Sq: 

{So<f>,  ^)  =  E  ( tiH\<l>?^r  +  H{<t^^  -  <f^-'^f)  , 
t=l   J=l 

where  4>\  \  (f)\  \  4>i  =  <f>\  are  the  values  of  the  function  <^''  G  Vo  at  the  nodes  of 
the  triangle  T,.  The  dimension  of  the  operator  So  is  equal  to  the  number  of  the 
nodes  which  belong  to  S^  \  Tq.  Let  tq  be  a  piecewise  linear  interpolation  operator: 

To:  Wo  ^  W„2  . 

According  to  Lemma  2.2,  we  can  set 

Eq     =  To5o     Tq     . 

Let  us  now  consider  the  subspaces  T^,  i  =  1, . . . ,  k.  From  Lemmas  2.1,  2.2  and  the 
definition  of  V,,  we  can  define 

where  S,,;v,  5",,^  generate  norms  in  Wi/2,i  and  W  1/2,,,  respectively.  To  define  5,,Ar 
and  S'^-p  we  will  again  use  Lemma  2.3.  Let  T,  be  embedded  into  the  square  IT,  and 
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the  triangulation  of  T,  be  complemented  to  triangulation  U'^  of  the  larger  region. 
In  the  square  Hi,  we  will  also  consider  the  auxiliary  square  mesh  Q'^  with  the  step 
ho:  ceh  <  ho  <  /imin/\/2-  Let  W{Q'-)  be  a  space  of  mesh  functions  which  is  defined 
on  Q'^  and  let  Sq,  be  an  operator  which  generates  an  equivalent  norm  in  the  trace 
space  W{Q'^).  According  to  Lemma  2.3,  we  can  define 

where 

rn,u,  =  p,ii,o,)p: . 

Here  /,  is  the  identity  matrix  which  corresponds  to  nodes  belonging  to  T,,  O,  the 
nuUmatrix  which  corresponds  to  nodes  belonging  to  IT,  \  T,,  and  P;  a  permutation 
matrix.  The  operators  Rv,n,,  Rq,  are  defined  exactly  as  were  those  used  when 
constructing  Bq,,  in  the  case  of  three  dimensions.  The  operator  5i,p  is  the  extension 
of  the  operator  S,,-p  by  zero.  We  can  thus  define 

E-^    =  E+  +  E+  +  •  •  •  +  E+  , 

Bt  =  tHi:-X  ,  (4-1) 

The  following  theorem  is  valid. 

Theorem.   There  exist  positive  constants  c-j,  cg,  independent  of  h  and  H ,  such 
that 

c-r{Bu^u)  6  (Au,ti)  G  Cs{Bu,u)  . 

Here  A  is  defined  by  (1.3)  and  B  is  defined  by  (4.1). 
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